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Abstract We review a recently developed method, based on a pseudoparticle representation 
of correlated electrons, to describe both Fermi liquid and non-Fermi liquid be- 
havior in quantum impurity systems. The role of the projection onto the physical 
Hilbert space and the impossibility of slave boson condensation are discussed. 
By summing the leading coherent spin and charge fluctuation processes in a fully 
self-consistent and gauge invariant way one obtains the correct infrared behavior 
of the pseudoparticles. The temperature dependence of the spin susceptibility 
for the single channel and two-channel Anderson models is calculated and found 
to agree well with exact results. 



1. INTRODUCTION 

Highly correlated electron systems are characterized by a strong repulsion 
between electrons on the same lattice site, effectively restricting the dynamics to 
the Fock subspace of states without double occupancy of sites. The prototype 
model for such systems is the Anderson impurity model, which consists of 
an electron in a localized level Ed < (called d-level in the following) with 
on-site repulsion U, hybridizing via a transition matrix element V with one 
or several degenerate conduction electron bands or channels [1]. Depending 
on the number of channels M, the model exhibits the single- or the multi- 
channel Kondo effect, where at temperatures T below the Kondo temperature 
Tk the local electron spin is screened (M = 1) or overscreened (M > 2) by the 
conduction electrons, leading to Fermi liquid (FL) or to non-Fermi liquid (NFL) 
behavior [2] with characteristic low-temperature singularities, respectively. 
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As perhaps the simplest model to investigate the salient features of correla- 
tions induced by short-range repulsion, the Anderson model plays a central role 
for the description of strongly correlated electron systems: In the limit of large 
spatial dimensions [3] strongly correlated lattice systems reduce in general to 
a single Anderson impurity hybridizing with a continuum of conduction elec- 
tron states whose properties are determined from a self-consistency condition 
imposed by the translational invariance of the system [4]. Quantum impurity 
models have received further interest due to their relevance for mesoscopic 
systems like single electron transistors or defects in quantum point contacts. 

The above-mentioned systems call for the development of accurate and flex- 
ible theoretical methods, applicable to situations where exact solution methods 
are not available. We here present a general, well-controlled auxiliary boson 
technique which correctly describes the FL as well as the NFL case of the 
generalized SU(N)xSU(M) Anderson impurity model. As a standard diagram 
technique it has the potential to be generalized for correlated lattice problems 
as well as for non-equilibrium situations in mesoscopic systems. 

2. AUXILIARY PARTICLE REPRESENTATION 
2.1 SU(N) x SU(M) ANDERSON IMPURITY MODEL 

The auxiliary particle method [5] is a powerful tool to implement the effective 
restriction to the sector of Fock space with no double occupancy imposed by a 
large on-site repulsion U. The creation operator for an electron with spin a in 
the d-level is written in terms of fermionic operators f a and bosonic operators b 
as d' a = f^b. This representation is exact, if the constraint that the total number 
operator of auxiliary fermions f a and bosons b is equal to unity is obeyed. /J 
and ftt may be envisaged as creating the three allowed states of the impurity: 
singly occupied with spin a or empty. 

In view of the possibility of both FL and NFL behavior in quantum impurity 
systems mentioned in the introduction it is useful to introduce M degenerate 
channels for the conduction electron operators cj.„, labeled fj, = 1,2, ... , M, 
in such a way that in the limit of impurity occupation number — > 1 (Kondo 
limit) the M-channel Kondo model is recovered, i.e. the model obeys an 
SU(M) channel symmetry. The slave bosons then form an SU(M) multiplet bp_ 
which transforms according to the conjugate representation of SU(M), so that 
\i is a conserved quantum number. Generalizing, in addition, to arbitrary spin 
degeneracy N, a = 1,2,...,N, one obtains the SU(N) x SU(M) Anderson 
impurity model in pseudoparticle representation 



H = J2 £ K> C *> + E dT,tif° + V J2 (%M + ' (L1) 
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where the local operator constraint Q = J2a ftfa- + bpbp, = 1 must be 
fulfilled at all times. 

2.2 GAUGE SYMMETRY AND PROJECTION ONTO 
THE PHYSICAL FOCK SPACE 

The system described by the auxiliary particle Hamiltonian (1.1) is invariant 
under simultaneous, local U(l) gauge transformations, f a — > fcr^ T \ bp — > 
bp,e l ^ T \ with <p(r) an arbitrary, time dependent phase. While the gauge 
symmetry guarantees the conservation of the local, integer charge Q, it does 
not single out any particular Q, like Q = 1. In order to effect the projection 
onto the Q = 1 sector of Fock space, one may use the following procedure 
[6, 7]: Consider first the grand-canonical ensemble with respect to Q and the 
associated chemical potential —A. The expectation value in the Q = 1 subspace 
of any physical operator A acting on the impurity states is then obtained as 

{A) = lim _ lim #l , a2) 



A-,00 9_tr[e-' 3 ( H + x Q)} G a^oo (Q) 



G 



where the index G denotes the grand canonical ensemble and ( is the fugacity 
£ = e~ l3X . In the second equality of Eq. (1.2) we have used the fact that 
any physical operator A acting on the impurity is composed of the impurity 
electron operators d a , dj., and thus annihilates the states in the Q = sector, 
A\Q = 0} = 0. It is obvious that the grand-canonical expectation value in- 
volved in Eq. (1.2) may be factorized into auxiliary particle propagators using 
Wick's theorem, thus allowing for the application of standard diagrammatic 
techniques. The local U(l) gauge symmetry, which must not be broken ac- 
cording to Elitzur's theorem [8], precludes the existence of an auxiliary Bose 
condensate. The latter would imply a spurious phase transition at finite T in 
quantum impurity models and should, therefore, be avoided in any approxima- 
tion. 

It is important to note that, in general, A plays the role of a time dependent 
gauge field. In Eq. (1.2) a time independent gauge for A has been chosen. 
In this way, the projection is only performed at one instant of time, explicitly 
exploiting the conservation of the local charge Q. Thus, choosing the time 
independent gauge means that in the subsequent development of the theory, the 
Q conservation must be implemented exactly. This is achieved in a systematic 
way by means of conserving approximations [9], i.e. by deriving all self- 
energies and vertices by functional derivation from one common Luttinger- 
Ward functional <I> of the fully renormalized Green's functions, 



£&,/, c = S^{G b ,G f ,G c }/SG b j tC . 



(1.3) 
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This amounts to calculating all quantities of the theory in a self-consistent way, 
but has the great advantage that gauge field fluctuations need not be considered. 

2.3 INFRARED THRESHOLD BEHAVIOR OF 
AUXILARY PROPAGATORS 

The projection onto the physical subspace, Eq. (1.2), implies that the pseudo- 
fermion and slave boson Green's functions Gf, G b are definied as the usual 
time-ordered, grand canonical expectation values of a pair of creation and 
annihilation operators, however evaluated in the limit A — > oo. It follows that 
the traces involved in Gf, G5 are taken purely over the the Q = sector 
of Fock space, and thus the backward-in-time contribution to the auxiliary 
particle propagators vanishes. Consequently, the auxiliary particle propagators 
are formally identical to the core hole propagators appearing in the well-known 
X-ray problem [11], and the long-time behavior of Gf (G5) is determined 
by the orthogonality catastrophe [12] of the overlap of the Fermi sea without 
impurity (Q = 0) and the fully interacting conduction electron sea in the 
presence of a pseudofermion (slave boson) (Q = 1). It may be shown that 
the auxiliary particle spectral functions have threshold behavior with vanishing 
spectral weight at T = for energies uj below a threshold E Q , and power law 
behavior above E a , Af^(u) oc 0(w — E )u~ a f' b . 

For the single-channel Anderson model, which is known to have a FL ground 
state, the threshold exponents may be deduced from an analysis in terms of 
scattering phase shifts, using the Friedel sum rule, since in the spin screened 
FL state the impurity acts as a pure potential scatterer [13, 14, 15, 16], 



a/= 2n^_n| ; °6 = 1-^ (N>1,M = 1) (1.4) 

These results have been confirmed by numerical renormalization group (NRG) 
calculations [17] and by use of the Bethe ansatz solution in connection with 
boundary conformal field theory (CFT) [18]. On the contrary, in the NFL case 
of the multi-channel Kondo model the threshold exponents have been deduced 
by a CFT solution [19] as 

M N 
at = , a b = (N > 2, M > N) (1.5) 

Since the dependence of a/, a b on the impurity occupation number n<i shown 
above originates from pure potential scattering, it is characteristic for the FL 
case. The auxiliary particle threshold exponents are, therefore, indicators for 
FL or NFL behavior in quantum impurity models of the Anderson type. 
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Figure 1.1 Auxiliary fermion and boson spectral functions Af(uj), Ab(u>) near uj = 0, evalu- 
ated within NCA in A saddlepoint approximation for T = W~ 6 D, e d = -0.8D, F = 0.15D, 
Tk — 7 ■ 10~ 5 D; D = half band width. Power law behavior for uj — > and spectral weight at 
ui < are seen. The inset shows At(uj) at larger uj. 

2.4 SADDLEPOINT PROJECTION 

Equivalently to the method shown in Eq. (1.2) the projection onto the Q = 1 
subspace can be represented as an integration over the parameter A, 

(A) = J V[f, /; b, b; c, c] J ^ e-/^(Q-iHW;M;c,c-}] {A)q (L6) 

where S is the action of the Anderson model corresponding to Eq. (1.1). This 
representation opens up the possibility of constructing an approximate projec- 
tion scheme [10]: Performing the A-integration in saddlepoint approximation 
(ASPA), i.e. evaluating Eq. (1.6) at the stationary point of the grand canonical 
action So = S + i\(Q — 1), is equivalent to determining — iX as a (real) 
thermodynamical chemical potential A„, fixing (Q) = 1. The ASPA respects 
several of the exact properties discussed in sections 2.2 and 2.3: (1) Although 
the threshold property of the auxiliary spectral function no longer holds, their 
infrared powerlaw behavior is preserved. Within NCA (section 3.1) it can be 
shown explicitly that the ASPA does not change the exponents a/, a&, and the 
same may be expected to hold for higher order self-consistent approximations 
(section 3.3). (2) Since the boson field is treated as a purely fluctuating field, 
Bose condensation does not occur in the ASPA. Technically, this is achieved in 
that the boson spectral function Aj,(uj), not restricted by the threshold property 
in this approximation, acquires weight at negative freqencies, which is negative 
because of the stability requirement ujA^lo) > for bosonic functions. The 
latter is seen in Fig. 1.1. Obviously, the saddlepoint projection is immediately 
generalizable to lattice problems, since the average local charge (Q) and thus 
A are space independent. 




6 



o \i O-co a' |i Sl-a 

T (cf) 
X CO X ' CO' 

Figure 1.2 Diagrammatic representation of the Bethe-Salpeter equation defining the conduc- 
tion electron-pseudofermion T-matrix T^'. 

3. CONSERVING SLAVE PARTICLE T-MATRIX 
APPROXIMATION 

3.1 NON-CROSSING APPROXIMATION (NCA) 

The conserving formulation discussed in section 2.2 precludes mean field 
approximations which break the U(l) gauge symmetry, like slave boson mean 
field theory. Although the latter can in some cases successfully describe the 
low T behavior of models with a FL ground state, it leads to a spurious phase 
transition at finite T and, in particular, fails to describe NFL systems. 

Rather, the approximation should be generated from a Luttinger-Ward func- 
tional <I>. Using the hybridization V as a small parameter, one may generate 
successively more complex approximations. The lowest order conserving ap- 
proximation generated in this way is the Non-crossing Approximation (NCA) 
[20, 21], defined by the first diagram in Fig. 1.3, labeled "NCA". The NCA 
is successful in describing Anderson type models at temperatures above and 
around the Kondo temperature Tk, and even reproduces the threshold expo- 
nents Eq. (1.5) for the NFL case of the Anderson impurity model. However, 
it fails to describe the FL regime at low temperatures. This may be traced 
back to the failure to capture the spin-screened Kondo singlet ground state of 
the model, since coherent spin flip scattering is not included in NCA, as seen 
below. 

3.2 DOMINANT LOW-ENERGY CONTRIBUTIONS 

In order to eliminate the shortcomings of the NCA mentioned above, we may 
use as a guiding principle to look for contributions to the vertex functions which 
renormalize the auxiliary particle threshold exponents to their correct values, 
since this is a necessary condition for the description of FL and NFL behavior, as 
discussed in section 2.3. As shown by power counting arguments [22], there are 
no corrections to the NCA exponents in any finite order of perturbation theory. 
Thus, any renormalization of the NCA exponents must be due to singularities 
arising from an infinite resummation of terms. In general, the existence of 
collective excitations leads to a singular behavior of the corresponding two- 
particle vertex function. In view of the tendency of Kondo systems to form a 
collective spin singlet state, we expect a singularity in the spin singlet channel 
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Figure 1.3 Diagrammatic representation of the Luttinger-Ward functional generating the 
CTMA. The terms with the conduction electron lines running clockwise (labelled "spin fluc- 
tuations") generate T (c ^', while the terms with the conduction electron lines running counter- 
clockwise (labelled "charge fluctuations") generate T^ cb \ The two-loop diagram is excluded, 
because it is not a skeleton. 



of the pseudofermion-conduction electron vertex function. It is then natural to 
perform a partial resummation of those contributions which, at each order in 
the hybridization V, contain the maximum number of spin flip processes. This 
amounts to calculating the conduction electron-pseudofermion vertex function 
in the "ladder" or T-matrix approximation, , where the irreducible vertex 
is given by V 2 G\ ) . The Bethe-Salpeter equation for reads (Fig. 1.2), 

Tl c T ^,{iu n ,iJ n ,i£L n ) = + V 2 G bfi (iuj n + iJ n - itt n )5 aT >5 Ta > 

- V 2 Tj2G bfl (iu;n + i^-in n )x (1.7) 

where a, r, a', t' represent spin indices and \i a channel index. A similar 
integral equation holds for the charge fluctuation T-matrix T^; it is obtained 
from by interchanging f a <-> and c afl «-> cj.„. Inserting NCA Green's 
functions for the intermediate state propagators of Eq. (1.7), we find at low 
temperatures and in the Kondo regime (n^O.7) a pole in the singlet channel 
of at an energy of the c-f pair which scales with Tk, indicating the onset 
of spin singlet formation. 

3.3 SELF-CONSISTENT FORMULATION 

In order to find a gauge invariant approximation incorporation multiple spin 
flip processes as discussed in the previous section, we have to find the Luttinger- 
Ward functional, which by second functional derivation w.r.t. G c and Gf or 
Gb generates the T-matrices and T( cb \ respectively. This generating 

functional <I> is shown diagrammatically in Fig. 1.3. The auxiliary particle self- 
energies are obtained in the conserving scheme as the functional derivatives 
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Figure 1.4 The fermion and boson threshold exponents af, a b are shown for JV = 2, M = 1 
in dependence of the average impurity occupation nd- Solid lines: exact values, Eq. (1.4); 
Symbols with error bars: CTMA; dashed lines: NCA. 

of with respect to Gf or G&, respectively (Eq. (1.3)), and are, in turn, 
nonlinear functionals of the full, renormalized auxiliary particle propagators. 
This defines a set of self-consistency equations, which we term conserving 
T-matrix approximation (CTMA). 

3.4 RESULTS 

We have solved the CTMA equations numerically for a wide range of im- 
purity occupation numbers rid from the Kondo to the empty impurity regime 
both for the single-channel and for the two-channel Anderson model down to 
temperatures of the order of at least W~ 2 Tk- 

The solution of the CTMA equations forces the T-matrices to have vanish- 
ing spectral weight at negative COM frequencies $7. Indeed, the numerical 
evaluation shows that the poles of and are shifted to SI = by 
self-consistency, where they merge with the continuous spectral weight present 
for SI > 0, and thus renormalize the threshold exponents of the auxiliary spec- 
tral functions. For N = 2, M = 1 the threshold exponents af, a.\, extracted 
from the numerical solutions are shown in Fig. 1.4. In the Kondo limit of the 
multi-channel case (N > 2, M = 2, 4) the CTMA solutions are found not 
to alter the NCA values and reproduce the the correct threshold exponents, 
a f = M/(M + N), a b = N/(M + N). 

The good agreement of the CTMA exponents with their exact values over 
the complete range of for the single-channel model and in the Kondo 
regime of the multi-channel model may be taken as evidence that the T-matrix 
approximation correctly describes both the FL and the non-FL regimes of 
the SU(N)xSU(M) Anderson model (N=2, M=l,2,4). The static spin sus- 
ceptibility x of the single- and of the two-channel Anderson model in the 
Kondo regime calculated within CTMA as the derivative of the magnetization 
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Figure 1.5 Static susceptibility of the single-channel (N = 2, M = 1) and the two-channel 
(JV = 2, M = 2) Anderson impurity model in the Kondo regime (_B d = -0.8 A V = O.LD, 
Lande factor # = 2). In the single-channel case, the CTMA and NCA results are compared to 
the NRG result (T = 0, same parameters) [23] and to the Bethe ansatz for the Kondo model 
[24]. The CTMA susceptibility obeys scaling behavior in accordance with the exact results (not 
shown). 

M = i(//x^(n/| —nfi) with respect to a magnetic field H is shown in Fig. 1.5. 
Good quantitative agreement with exact solutions is found for N = 2, M = 1 
(FL). For N = 2, M = 2 (NFL) CTMA correctly reproduces the exact [25] 
logarithmic temperature dependence below the Kondo scale Tk- In contrast, 
the NCA solution recovers the logarithmic behavior only far below Tk- 



We have reviewed a novel technique to describe correlated quantum im- 
purity systems with strong onsite repulsion, which is based on a conserving 
formulation of the auxiliary boson method. The conserving scheme allows to 
implement the conservation of the local charge Q without taking into account 
time dependent fluctuations of the gauge field A. By including the leading 
infrared singular contributions (spin flip and charge fluctuation processes), 
physical quantities, like the magnetic susceptibility, are correctly described 
both in the Fermi and in the non-Fermi liquid regime, over the complete tem- 
perature range, including the crossover to the correlated many-body state at 
the lowest temperatures. As a standard diagram technique this method has the 
potential to be applicable to problems of correlated systems on a lattice as well 
as to mesoscopic systems out of equilibrium via the Keldysh technique. 

We wish to thank S. Bocker, T.A. Costi, K. Haule, S. Kirchner, A. Rosch, 
A. Ruckenstein and Th. Schauerte for stimulating discussions. This work is 
supported by DFG through SFB 195 and by the Hochleistungsrechenzentrum 
Stuttgart. 
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